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1.  INTRODUCTION  AND  SUMMARY 


1.1  INTRODUCTION 

Radomes  affect  the  performance  of  airborne  radar  and  guidance  systems  by  reducing 
guidance  accuracy  and  radar  range.  These  system  effects  limit  existing  missiles,  and 
they  restrict  the  development  of  newer  missiles  that  have  increased  speeds,  higher 
altitudes,  greater  frequency  bandwidths,  or  smaller  dimensions.  The  need  for  im- 
proved accuracy  requires  accurate  methods  for  radome  analysis.  Although  adequate 
methods  exist  for  designing  radomes  that  enclose  antennas  with  diameters  greater 
than  several  wavelengths,  these  methods  fail  for  smaller  antennas  and  radomes.  The 
development  of  small  or  broadband  radomes  is  largely  empirical  so  that  data  are 
available  only  after  a radome  is  fabricated  and  tested.  This  approach  is  tedious  and 
expensive,  and  it  delays  the  synthesis  of  systems.  A more  general,  analytical  design 
method  is  desirable. 

The  development  of  a general  design  method  requires  a model  that  relates  radome 
composition  and  configuration  to  system  effects.  The  link  is  the  physical  effects  of  a 
radome;  namely,  boresight  error,  attenuation,  and  sidelobe  level  increases.  These 
three  effects  can  be  qualitatively  and  quantitatively  derived  from  farfield  patterns  of 
the  antenna  in  the  radome.  * Consider  sum  mode  patterns  with  and  without  radome  as 
in  Figure  1-1.  Attenuation  or  loss  is  the  reduction  in  intensity  at  the  peak  of  the  main 
beam,  and  boresight  error  is  the  angular  shift  of  the  beam  peak.  For  monopulse  or 
conical  scan  antennas,  boresight  error  is  the  angular  shift  of  the  tracking  minimum 
Instead  of  the  peak.  Sidelobe  levels  changes  are  the  changes  in  the  intensity  of  sub- 
sidiary maxima  in  the  antenna  pattern. 

Antenna  patterns  are  computed  by  evaluating  diffraction  integrals.  We  assume 
receiving  operation  in  order  to  consider  plane  waves  incident  on  the  radome.  The 
complex-valued  amplitude  of  a rectangular  component  in  the  radome-bounded  region  is 

ER  = jj FTEJdA  , 

where  E*  is  the  incident  field  amplitude,  F is  the  receiving  aperture  distribution,  and 
T is  the  radome  transmittance,  which  is  a complex-valued  function  of  position  for  fixed 
frequency,  wave  polarization,  and  incidence  direction.  The  integration  extends  over 
the  region  in  w MF  has  nontrivial  magnitude.  The  receiving  power  is  assumed  pro- 
portional to  |E  | . 

♦The  pattern  depends  on  the  orientation  of  the  antenna  in  the  radome  as  well  as  fre- 
quency and  wave  polarization  direction. 
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The  transmittance  T In  Equation  (1)  describes  the  electromagnetic  effects  of  the 
radome.  If  we  write  T = T e1^*  , A$  is  the  phase  delay  produced  by  the  radome. 

The  quantity  A $ Is  significant  because  It  determines  equlphase  surfaces,  and  the  tilt 
of  these  surfaces  from  the  plane,  Incident  wave  Is,  to  first  order,  proportional  to 
boreslght  error.  Figure  1-2  suggests  a distorted  wavefront.  Because  T strongly 
Influences  boreslght  error.  Its  accurate  computation  Is  necessary  for  accurate  per- 
formance predictions. 

The  transmittance  depends  on  several  electromagnetic  wave  processes.  The  known 
mechanisms  are  suggested  In  Figure  1-3  for  externally  incident  plane  waves.  The 
direct  wave  usually  has  the  largest  magnitude.  For  large  radomes  and  antennas  the 
direct  wave  can  be  accurately  computed  over  most  of  the  aperture  area  by  approxi- 
mating the  radome  as  locally  plane.  However,  this  approximation  fails  in  the  shadow 
of  the  tip  because  the  radius  of  curvature  of  the  radome  tip  approximates  the  wave- 
length. Some  evidence  exists  for  a scattered  wave  centered  at  the  tip  or  the  edge  of 
a wedge.3 

Data  for  other  radomes  suggests  that  if  a wave  is  scattered  by  the  tip,  the  wave  is 
very  directional;  a simple  cylindrical  or  spherical  wave  may  be  an  overly  approximate 
assumption.  Any  tip  scattered  wave  becomes  more  significant  as  the  diameter  of  the 
receiving  antenna  is  reduced  because  the  projection  of  the  tip  region  becomes  a larger 
fraction  of  the  aperture  area.  Tip  diffraction  may  be  more  important  for  direction- 
finding antennas  that  sense  the  difference  between  fields  for  halves  of  the  aperture;  for 
larger  angles  of  arrival,  the  direct  waves  encounter  large  radii  of  curvature  and  pro- 
duce balanced  fields  that  cancel,  but  a tip  wave  would  have  distinct  phase  values. 

The  reflected  wave,  which  comes  from  the  shadowed  side  of  the  radome,  occurs  for 
incidence  directions  sufficiently  far  from  axial  incidence.  The  angle  is  approximately 
that  of  a cone  tangent  to  the  radome  because  the  incident  wave  must  first  propagate 
through  the  directly  illuminated  side.  The  reflected  wave  increases  sidelobe  levels, 
and  it  causes  boresight  error  for  larger  gimbal  angles. 

Guided  waves  exist  in  radomes.  The  modes  are  hardly  cataloged  for  any  special  case. 
In  general,  evidence  exists  for  slab  guided  modes  in  wedges  and  cones.  Near  the  tip, 
the  radome  may  be  a hollow  circular  waveguide.  Guided  modes  are  negligible  in  many 
applications  if  the  aperture  is  large,  but  they  are  growing  in  importance  for  newer 
smaller  seekers  that  have  antennas  with  small  clearance  from  the  radome.  In  addition, 
a guided  wave  can  excite  a tip  diffracted  wave  and  scattering  at  the  base. 
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Until  recently  the  direct  wave  was  the  main  concern.  However  the  other  wave  proces- 
ses, which  had  produced  rather  subtle  effects  on  systems,  have  become  more  signi- 
ficant as  radome  boresight  error  slope  requirements  become  more  stringent.  In 
addition  low  sidelobe  antennas  require  accurate  design  methods  especially  for  reason- 
able bandwidths. 

This  report  describes  analyses  and  experiments  on  wave  propagation  through  hollow 
dielectric  shells.  The  purpose  was  to  better  describe  the  wave  processes  shown  in 
Figure  1-3.  Our  aim  was  to  analyze  wavefront  structure  and  the  excitation  of  the 
various  wave  processes  to  estimate  magnitude  of  the  constituent  fields.  The  results 
have  helped  develop  anisotropic  structures  to  reduce  boresight  error  and  its  variation 
with  wave  polarization.  Although  completely  polarization  independent  radomes  seem 
unlikely,  some  reduction  of  this  problem  seems  possible  with  dielectric  gratings  or 
periodically  perforated  metallic  structures.  Finally  the  results  of  the  research  are 
expected  to  help  the  analysis  of  high  temperatures  on  radome  boresight  error. 
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Figure  1-1.  Radome  Effects  from  Antenna  Patterns 
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» INCIDENT  WAVEFRONT 


Figure  1-2.  Distorted  Wavefront 
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Figure  1-3.  Wave  Mechanisms  in  Hollow  Dielectric 


1.2  SUMMARY 

This  report  describes  analyses  and  experiments  on  wave  propagation  through  hollow 
dielectric  shells.  The  purpose  was  to  describe  quantitatively  the  wave  process  in 
radomes^  see  Figure  1-3.  -Specifically  we  sought  to  analyze  wavefront  structure  and  to 
estimate  the  magnitudes  of  the  elementary  waves  that  combine  to  form  the  wavefront. 
The  research  applies  to  the  design  of  new  missile  radomes,  to  the  analysis  of  aero- 
dynamic heating  on  radome  performance,  and  to  the  development  of  anisotropic  radomes 
that  reduce  the  polarization  dependence  of  radome  performance.  v 

Section  2 summarizes  the  theory,  which  is  the  moment  method.  It  is  restricted  to  two- 
dimensional  problems  in  its  present  form.  Section  3 describes  calculated  results  for 
plane  waves  incident  on  flat  slabs  and  on  hollow  wedges  for  a range  of  incidence  direc- 
tions. Emphasis  is  on  perpendicular  polarization,  but  data  for  parallel  also  are  given. 
Results  are  given  for  several  slab  and  wedge  dimensions.  Section  4 describes  near- 
field experiments.  The  results  are  compared  with  computed;  and  an  accuracy  criterion 
is  given.  Section  5 interprets  the  results.  The  measured  and  computed  results  show 
guided  waves  are  strongly  excited  on  slabs  and  wedges.  The  moment  method  correctly 
gives  the  magnitude  of  excitation.  An  independent  model  predicts  spacings  of  fringes 
formed  by  interfering  guided  and  refracted  waves.  Scattering  by  edges  also  occurs. 
These  processes  plus  propagation  through  and  reflection  from  the  flat  sides  of  a wedge 
produce  complicated  fields.  A connection  was  made  between  observed  boresight  error 
properties  of  axially  symmetric  radomes  and  guided  waves. 
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2. 


ANALYSIS 


2.1  INTRODUCTION 

Radome  analysis  requires  approximations  because  few  radomes  are  bounded  by  sur- 
faces that  are  complete  co-ordinate  surfaces  in  separable  systems.  Cylinders, 
spheres,  and  large  flat  sheets  can  be  analyzed  by  boundary  value  methods;  however, 
ogival  radomes  usually  are  analyzed  with  diffraction  integrals,  and  the  radome  is 
described  as  locally  flat.  Although  approximate  methods  are  often  adequate,  their 
accuracy  is  unknown  in  many  cases.  Discrepancies  between  measured  and  computed 
quantities  are  largest  for  small  antennas  and  small  radomes. 

The  causes  of  the  discrepancies  are  only  partially  known.  The  approximation  of  local 
planeness  in  computing  transmittance  is  mentioned  by  several  authors.  Our  experience 
suggests  additional  causes,  namely  undersampling  in  computing  transmittance.  To 
appreciate  undersampling  consider  Figure  2-1A.  A single  incident  wave  normal,  a 
direct  ray,  passes  through  a point  in  the  radome -bounded  region.  Transmittance  is 
computed  by  approximating  the  radome  as  locally  flat.  This  method  gave  good  accura- 
cy when  the  antenna  in  the  radome  had  diameter  10  wavelength,  but  it  failed  for  an 
antenna,  diameter  wavelength  in  a smaller  radome.  In  contrast,  the  method  suggested 
by  Figure  2- 1C  gave  excellent  accuracy  for  the  smaller  radome.  We  interpret  these 
results  to  suggest  that  the  radome  surface  must  be  adequately  sampled.  Sufficiently 
dense  sampling  introduces  the  variation  of  the  surface  normal  into  a sum  of  terms  that 
describe  transmittance  at  a point.  The  lack  of  this  variation  in  the  direct  ray  method 
seems  a cause  of  discrepancies. 

Another  cause  of  discrepancies  seems  to  be  the  omission  of  some  wave  processes  from 
most  analytical  methods.  As  mentioned  in  Section  1,  guided  waves,  scattering  by  the 
tip,  and  by  the  base  of  the  radome  exist,  but  these  waves  are  usually  omitted  from 
numerical  analyses.  A main  objective  of  this  project  was  to  determine  the  magnitude 
of  the  distinct  wave  types  in  radome-bounded  regions. 


2-1 


Because  accurate  radome  analysis  Is  involved,  we  considered  simpler  problems,  the 
hollow  dielectric  wedge  and  the  flat  slab,  both  with  finite  width,  as  shown  in  Figure  2-2. 
In  addition,  we  considered  the  half  cylinder.  These  problems  are  preliminaries  to 
axially- symmetric  radomes.  We  considered  two  directions  of  wave  polarization.  One 
is  TM  or  perpendicular  polarization.  The  other  is  TE  or  parallel  polarization.  The 
approach  is  to  utilize  the  moment  method  for  the  numerical  solution  of  integral 
equations. 

The  theoretical  basis  for  the  analysis  of  slabs,  hollow  wedges,  and  hollo^  flinders  is 
the  moment  method;  in  particular,  a form  developed  by  J.  H.  Richmond.  * The  next 
two  subsections  summarize  the  theory. 

2.2  TM  CASE,  OR  PERPENDICULAR  POLARIZATION 

The  starting  point  is  an  integral  equation  for  the  scattered  electric  field.  This  equation 
is  derived  from  Maxwell’s  equations  with  the  definition  ti^at  E , the  total  field,  which 
is  the  field  in  the  presence  of  a scatterer^Ls  the  sum  of  E , the  incident  field,  which  is 
the  field  without  the  scatterer,  and  the  E the  scattered  field.  That  is 

T I S 

E = E + E . (1) 

Maxwell’s  equations  are  then  rewritten  in  the  form  of  the  free  space  equations  with  a 
polarization  current  that  describes  the  effect  of  the  dielectric  scatterer.  For  perpen- 
dicular polarization  the  integral  equation  gives  the  scattered  field  as 

ES(x,y)  = (ik2/4)  //  ( /c-l)ET(x',y')Ho(1)  (kp  )dx’dy',  (2) 

2 2 1/2 

where  P = [(x-x’)  + (y-y')  ] , with  (x,y)  the  co-ordinates  of  the  observation 

point  and  (x',y')  those  of  the  integration  point.  H ' '(k  P)  is  a zero  order  Handel 
function  of  the  first  kind,  it  represents  an  outwarc?  going  wave  witi^our  assumed  time 
convention  exp  (-iw  t),  with  t time  and  frequency  a:.  By  adding  E in  Equation  (2)  to 
E1,  we  obtain  an  integral  equation  for  ET. 

T 2 T ill  I 

E - (ik  /4)  //(  k-1)E  (x',y')HoV  ; (kP)dx’dy'  = E . (3) 

To  solve  Equation  (3),  divide  the  two  dimensional  dielectric  scatterer  into  oplls  like 
those  in  Figure  2-3.  The  cells  must  be  small  enough  to  justify  assuming  E constant 
in  any  one  cell.  We  assume  the  relative  dielectric  constant  k is  constant  in  any  one 
cell.  From  Equation  (3)  ve  get 

T 2 m I 

E - (ik74)  //  ( « -l)E(x  ,y  )Hl  '(kP)dx'dy'  = E . (4) 

m m mm 
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In  Equation  (4)  m is  an  index  that  labels  the  cells,  and  N is  the  total  number  of  cells. 

In  fact  Equation  (4)  is  a system  of  N equations.  The  next  step  is  to  evaluate  the  double 
integrals  in  Equation  (4).  This  integration  is  straightforward  for  cells  that  are  right 
circular  cylinders  of  radius  a and  with  Infinite  length  (in  the  z*  direction,  which  Is 
perpendicular  to  the  x',y'  plane).  There  are  two  cases  to  consider.  In  one,  the  ob- 
servation point  is  outside  the  cell.  If  n denotes  the  integration  cell  and  m the  obser- 
vation cell,  then  for  m * n 

2 m 

(-ik  / 4)1  = (-i7rka/2)Ji(ka)HQV  (kpmQ),  (5) 


where  p denotes  the  distance  between  (x  , y ) and  xq,  y ) , I is  the  double  integral 
in  Equationn(4),  and  J denotes  the  first  order  Vessel  function.  In  the  second  case, 
m = n 


(-ik  /4)  I = (-1/2) 
With  Equations  (5)  and  (6),  we  obtain 


j^rk  a Hx(1)(ka)  + 2i  . 


N 


I 

n =1 


C E 
mn 


T 

n 


f 


(6) 


(7) 


for  1 s.  m«  N.  The  C are  matrix  elements. 

mn 

To  find  E it  is  necessary  to  solve  the  system  of  Equation  (7).  This  solution  may  be 
found  by  Inverting  the  matrix  (C  ).  The  matrix  inversion  is  restricted  by  storage 
capacity  of  computers  because  tfien assumption  that  E is  constant  in  any  one  cell 
requires  small  cells.  The  small  dimensions  require  large  N for  structures  compar- 
able in  size  to  practical  radomes. 

Direction  solution  of  equation  7 is  being  studied.  The  next  section  gives  calculations 
for  four  special  cases  that  were  chosen  to  reveal  some  aspects  of  wave  propagation  in 
and  around  radomes. 

2.3  TE  CASE  OR  PARALLEL  POLARIZATION 

This  case  is  more  complicated  than  the  preceding  because  another  term  arises  in  the 
Integral  equation  since  the  gradient  of  the  scalar  potential  contributes  to  the  electric 
field, 

9A 

E = - VP-  — • (8) 

o l 
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With  the  co-ordinate  system  of  Figure  2-2,  the  incident  field  is 


E1  = x E1  (x, y)  + yE1  (x,y); 
x y 


(9) 


The  total  field  E has  a similar  representation.  Two  sets  of  coupled  equations  result, 
we  obtain 


N 

2 


n=l 


(A  E + B E ) 
mn  xn  mn  yn 


xm 


N 


2 


(C  E + D E ) 
mn  xn  mn  yn 


(10) 


where  A , 
_ . mn, 
Reference  4. 


B , 


mn 


mn 


D are  matrix  elements;  their  definitions  are  given  in 
mn 
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PLANE  PARALLEL  TO 
INCIDENT  WAVEFRONT 
AS  INTEGRATION  SURFACE 


Figure  2-1.  Three  Methods  for  Computing  Radome  Boreslght  Error 


Figure  2-2.  Orientations  of  Slabs  and  Hollow  Wedge. 
The  dashed  lines  show  where  the  field  was  evaluated. 


For  the  wedge,  a range  of  incidence  directions  was  considered. 
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3. 


COMPUTATION 


This  section  gives  the  results  of  moment  calculations  for  slabs,  hollow  wedges,  and 
hollow  hemi-cyllnders.  The  co-ordinate  systems  of  Figures  3-1,  3-2,  and  3-3  were 
utilized.  Cell  radius  a was  chosen  to  make  the  sum  of  the  areas  of  the  approximating 
cylinders  equal  that  of  the  scattering  object.  Frequency  was  9.  375  GHz.  Dielectric 
constant  k was  2. 6;  to  simulate  acrylic  plastics  which  we  utilized  in  experiments. 
Loss  tangent  was  zero;  this  approximation  Is  good  because  the  value  for  acrylic 
plastics  is  0.  006. 

3. 1 PERPENDICULAR  POLARIZATION 

3. 1. 1 SLAB:  4 IN.  WIDTH,  NORMAL  INCIDENCE 

The  fields  inside  the  slabs  were  computed  at  the  center  of  the  approximating  cylinders 
to  study  numerical  convergence  and  to  compare  the  results  with  those  for  infinitely 
broad  sheets.  Incidence  was  normal  as  in  Figure  2-2a.  Figure  3-4  shows  I E I 
inside  a slab,  width  4 in.  and  thickness  0.  25  in.  Cell  radius  a was  0. 143  in. ; spacing 
was  0.  25  in.  The  number  of  cells  N was  17.  Note  the  valuq^  fluctuate  about  that  for 
infinite  flat  sheets.  For  comparison.  Figure  3-5  shows  lE  I inside  the  same  slab, 
but  for  N = 120  and  for  smaller  cells.  For  the  slab  midplane  the  results  are  similar 
for  both  values  of  N,  but  the  oscillations  have  greater  amplitude  in  Figure  3-5.  The 
field  also  depends  on  Y within  the  slab. 

T 12 

The  power  transmittance  (the  normalized  intensity  I E /E  I outside  the  slab)  of  the 
0.  25  in.  thick  slab  is  shown  in  Figure  3-6  for  N = 17  and  in  Figure  3-7  for  N = 120. 
These  figures  illustrate  how  transmittance  computations  depend  on  N. 

i T 12 

Figure  3-8  shows  I E /E  I for  a slab  with  thickness  0. 125  in.  In  this  case  the  slab 
was  approximated  with  a single  row  of  33  cylinders. 

3.  1.2  SLAB:  10  IN.  WIDTH,  NORMAL  INCIDENCE 

As  the  width  of  a slab  increases  we  expect  the  field  to  approach  that  computed  for  an 
infinitely  broad  slab.  Therefore  we  considered  a second,  broader  slab,  width  10  in. 
The  internal  field  magnitude  is  shown  in  Figure  3-9  for  two  thicknesses,  0. 125  and 
0.  250  in.  Figure  3-10  shows  transmittance  at  several  distances  behind  the  0.  25  in. 
slab;  all  are  in  the  nearfield  of  the  slab.  Figure  3-11  shows  transmittance  for  the 
0. 125  in.  slab.  Note  in  Figure  3-11  the  transmittance  differs  slightly  from  that  for 
a flat  slab  in  the  central  half  of  the  slab,  but  nearer  the  edges  the  moment  method 
result  oscillates  about  the  flat  sheet  value. 
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3.  1.  3 SLAB  WITH  DIELECTRIC  STRIPS 


Figure  3-12  shows  an  array  of  dielectric  cylinders  that  represents  thin  strips  on  a 
dielectric  slab.  Slab  thickness  was  0. 125  in. , and  width  was  4.  38  in.  Each  strip  was 
0.30  in.  high  and  0.060  in.  thick.  The  slab  was  represented  by  36  cylinders  each  with 
diameter  0.  071  in.  ^ and  each  strip  had  5 cylinders,  diameter  0.  034  in.  Computed 
values  of  I E/E1 1 “ are  shown  in  Figure  3-13  for  normal  incidence  on  the  slab  with 
strips. 

The  effects  of  the  strips  on  power  transmittance  is  illustrated  in  Figure  3-14,  which 
shows  the  computed  data  for  the  slab  with  and  without  the  15  strips.  The  strips  in- 
crease transmittance  averaged  over  the  slab  by  approximately  1/3  dB. 

3.1.4  SLAB,  45°  INCIDENCE 

For  45°  incidence,  as  in  Figure  2-2b,  the  internal  field,  at  the  centers  of  the  cylinders 
approximating  a slab,  is  shown  in  Figure  3-15  for  six  values  of  N(10,  11,  21,  41  and 
42)  corresponding  to  six  distinct  slab  lengths.  The  intensity  outside  the  slab  is  shown 
in  Figure  3-16  for  N = 41;  in  this  case  intensity  was  evaluated  on  the  line  0. 225  in  from 
the  slab  midplane  and  on  the  shadowed  side.  This  line  is  0.1  in.  from  the  slab.  The 
intensity  on  a transverse  line  parallel  to  the  incident  wavefront  is  shown  in  Figure  3-17 
for  two  values  of  N(41  and  42). 

The  internal  fields,  Figure  3-15  depend  on  N,  but  this  dependence  has  a physical 
rather  than  numerical  cause.  The  deep  minima  result  from  guided  waves  interfering 
with  internal  plane  wave  fields;  we  justify  this  interpretation  in  Section  5.  The  small 
rapid  oscillations  result  from  interference  of  guided  waves  travelling  in  opposite 
directions;  see  Section  5.  Addition  of  another  cylinder  increased  the  path  length  of 
one  of  the  waves  by  approximately  tr,  causing  maxima  and  minima  of  the  small 
oscillations  to  alternate  in  locations. 

Another  physical  cause  exists  for  fluctuations  on  transverse  paths  as  in  Figure  3-17. 
The  increase  of  N extends  the  slab.  For  small  N the  end  of  the  slab  is  closer  to  the 
source  than  the  observation  points.  The  larger  N the  end  is  farther  from  the  source. 
Figure  3-18  illustrates  this  point  through  calculations  of  the  scattered  Field. 

3. 1.  5 SLAB,  GRAZING  INCIDENCE 

For  grazing  incidence  as  in  Figure  2- 2c,  Figure  3-19  shows  the  intensity  of  the  total 
field  in  a 0. 25  in.  thick  10  in.  long  slab.  The  field  was  evaluated  at  the  ceop^s  of 
the  .cells  for  two  values  of  N(40  and  41).  We  see  from  the  4 dB  value  of  |E  I that 
| E | exceeds  the  unit  magnitude  of  the  incident  field. 


F 


The  external  Intensity^  jyi  Figure  3-20  for  two  lines  parallel  to  the  slab  plane,  and 
Figure  3-21  shows  I E I at  a larger  distance. 

3. 1.  6 WEDGE,  0.  25  IN.  THICK,  AXIAL  INCIDENCE 

T 

For  a hollow  wedge,  as  in  Figure  3-22  shows  Ie  I inside  the  walls  for  several  values 
of  N,  where  (N-l)/2  cells  are  in  each  wall  and  one  cell  is  at  the  tip.  The  field  on 
transverse  paths  (fixed  y values)  is  shown  in  Figure  3-23  for  the  wedge  represented 
by  two  slabs  of  equal  thickness,  and  in  Figure  3-24  for  two  slabs  of  slightly  different 
thickness.  Note  the  additional  variation  in  the  second  case.  The  field  variation  on 
the  symmetry  plane  is  shown  in  Figure  3-25. 


To  obtain  a comprehensive  view  of  the  field  we  drew  contour  plots  of  the  scatterec^ 
field;  see  Figure  3-26  for  magnitude  an<J  3-27  for  phase.  The  biggest  value  of  I E | 
was  1. 65  compared  to  the  unit  value  I E I . 


3. 1.  7 WEDGE,  0. 125  IN.  THICK  WALLS 

Figure  3-28a  and  b shows  the  field  intensity  inside  a hollow  90°  wedge,  thickness 
0. 125  in. , with  walls  10  in.  long;  results  are  given  for  four  incidence  directions  0 . 
equal  to  0°,  22.5°,  45°,  and  67.5°.  Note  0j  equal  0°  gives  axial  incidence:  0} 
equal  45°  gives  normal  incidence  on  one  slab  and  grazing  incidence  on  the  other.  The 
value  67.5°  gives  reflections. 

Because  many  graphs  are  necessary  to  represent  the  field  behind  a wedge,  we  pictorial 
ly  displayed  the  data.  Figures  3-29  and  3-30  show  these  plots  of  the  scattered  field 
E®  for  the  wedge  with  0. 125  in.  thick  walls.  The  data  were  quantized  into  11  levels 
so  the  representation  is  rather  coarse.  Therefore  we  increased  the  number  of  levels 
to  22;  see  Figure  3-31  and  3-32. 

3. 1.  8 HEMICYLENDER 


Figure  3-33  shows  the  geometry.  The  hemicylinder  had  diameter  4.  8 in.  and  thickness 
0.  125  in.  The  direction  of  the  incident  wave  was  specified  by  the  angle  0 . The  hemi- 
cylinder was  approximated  by  a set  of  small  solid  cylinders,  radius  a,  on  a hemicircle 
of  radius  R.  See  Figure  3-34.  Computations  were  done  for  R equal  2.  5 in.  and  a 
equal  0.  071  in. ; P the  distance  between  hemicylinders  was  . 125  in. 

T X 2 

Figure  3-35  shows  intensity  t E /E  I for  0 = 90°  at  two  values  of  y.  For  compari- 
son, the  figure  also  shows  values  computed  with  a ray  trgpin^  method  that  approxi- 
mates the  radome  as  locally  flat.  The  phase  delay  arg  E /E  is  shown  in  Figure  3-36. 
Figures  3-37  and  3-38  show  intensity  and  phase  for  = 45°. 
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Because  anisotropic  gratings  have  been  utilized  to  Improve  radomes,  we  applied  the 
moment  method  to  their  analysis.  Figure  3-39  shows  the  geometry.  Transmittance 
is  shown  in  Figure  3-40  for  o.  =0°. 

3.2  PARALLEL  POLARIZATION 

3.  2. 1 SLAB:  0. 125  IN.  THICK,  GRAZING  INCIDENCE 
T 

Figure  3-41  |E  I Insj^le  | s^ab,  thickness  0. 125  in.  length  4 in.  for  grazing  incidence. 
Figure  3-42  shows  IE  /E  I outside  the  slab  at  five  lines  parallel  to  the  slab. 

3.2.2  SLAB:  0.235  in.  THICK,  GRAZING  INCIDENCE 

T 

Figure  3-43  shows  lE  I inside  a slab,  thickness  0.235  in. , length  4 in.  for  grazing 
incidence.  Results  are  shown  for  ^ slijgl^  row  of  cylinders,  N = 17,  and  for  two  rows 
N = 68.  The  external  intensity  |E  /E  I “ is  shown  in  Figure  3-44. 
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Figure  3-4.  Internal  Field  Magnitude  for  1/4  in.  x 4 in.  x 18  in.  Slab,  k = 2.  6, 
for  Normal  Incidence.  Calculated  for  N = 17,  a = 0.  143  in.  , P = 0.  25  in.  , 
at  y = 0.  The  arrow  shows  the  value  predicted  by  flat  sheet  theory. 
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Figure  3-5.  Internal  Field  Magnitude  for  1/4  in.  x 4 in.  x 18  in.  Slab,  k = 2.  G,  for 
Normal  Incidence.  Calculated  for  N = 120,  a = 0.  053  in.  , p = 0.  100  in.  ; 

Hexagonal  Arrangement  of  Cells  for  y = 0 ( );  for  y = ±0.  08  in.  ( ).  The 

arrow  shows  the  value  predicted  by  the  theory  of  flat  sheets 
at  the  central  row  of  cylinders. 
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Figure  3-<5 . Power  transmittance  of  0.  24  b>  4.  0 in.  slab;  k = 2.  6;  0j  = 0*; 
Calculated:  N = 17,  a = 0.  143  in.  , P = 250  in 
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Figure  3-7.  Power  transmittance  of  slab  in  Figure  3-<i  'out  calculated  N ■ 120, 

a = 0. 053  in. , p = 0.  10  in. 


0 


E005 


2 h 


X uNCH' 


MEASURED  With  DiPOLE ‘ 

CALCULATED  ( I 


Figure  3-8.  Power  transmittance  of  0.125  by  4.0  in.  slab:  k = 2.C,  bj  = 908. 
N = 33,  a = 0.0715  in.,  P=  0.125  in. 


FL006 


0 1 1 1 1 1 L 

-4  -2  0 2 4 

X (INCH) 


Figure  3-9.  Total  Field  Amplitude  in  Dielectric  Slab,  k = 2.6,  Normal  Incidence. 
Computed  for  N=  41,  a = 0. 143  in. , P = 0.25  in.  (for  1/4  in.  thickness  the  dashed 
curve)  and  N=  81,  a=  0.0715  in. , P = 0.125  in.  (for  1/8  in.  thickness  the  solid  curve). 
The  arrows  indicate  the  values  predicted  by  flat  sheet  theory. 
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Figure  3-11.  Power  transmittance  of  0.125  by  10  by  18  in.  Plexiglas  slab; 

Calculated:  N = 81,  a=  .0715  in.,  P = . 125  in.  The  broken 
heavy  line  is  the  value  for  an  infinite  slab. 
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Figure  3-12.  Array  of  11  cylinders  to  represent  4.38  in.  wide  slab 
0.125  in.  thick  with  15  strips. 
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Figure  3-13.  Power  transmittance  of  0. 125  by  4.  38  in.  slab  with  15  strips:  Computed 
for  N = 111,  a = 0.  074  in.  and  0.  034  in.  , p = 0. 125  in.  and  0.  06  in. 
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Figure  3-14.  Power  transmittance  of  0. 125  by  4.  38  in.  slabs:  « = 2.  6,  flj  = 0°; 

with  15  dielectric  strips:  N = 111;  a = 0.  071  in. : 0.  034  in. 

p = 0. 125  in.  : 0.  060  in. 

without  strips:  N = 36:  a = 0.071  in.  ; p = 0.125  in. 


Figure  3-15.  Total  Field  Magnitude  Inside  Dielectric  Slabs  for  45°,  X = 1.259  in. 
Computed  for  a = 0.143  in.,  p = 0.250  in.  The  arrows  show 
theoretical  fringe  spacings  for  guided  waves. 
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Figure  3-16.  Intensity  in  Total  Field  at  Distance  0.  225  in.  from  Slab  Midplane. 
Thickness  0.  25  in. , length  10  in.  Computed  for  N = 41,  a = 0. 143  in. , p = 0.  25  in. , 
X = 1.259  in.,  k = 2.6.  The  numbers  between  minima  are  spacings  in  inches. 
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Figure  3-17.  Intensity  Behind  Dielectric  Slab,  k = 2.  6,  X = 1.  259  in.  for  45°  Incidence. 

Measured  for  0.  23  in.  x 10  in.  x IS  in.  slab  with  dipole  probe  at  x = 4. 48  in.  ( ). 

Calculated  for  N = 41  (°);  for  N = 42  (o).  P = 0.  25  in. , a = 0. 143  in. 
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Figure  3-19.  Total  Field  Intensity  Inside  Slab  1/4  in.  Thick,  10  in.  Long. 
N = 40  (— ),  N = 41  ( ),  a = 0. 143  in. , p = 0.  25  in. , k = 2.  6. 
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Figure  3-20.  Intensity  Near  Dielectric  Slab,  k = 2.  6,  1/4  in.  x 10  in.  x 18  in. , 
Calculated  for  N = 40,  a = 0. 143  in. , P = 0.  25  in. , \ = 1.  259  in. 
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Figure  3-23.  Intensity  Behind  Hollow  Dielectric  Wedge,  k = 2.  6,  With 
10  In.  x 18  In.  Walls,  Thicknesses  0.226  In.  and  0.232  In.  at  \=  1.259  In. 
Measured  with  dipole  probe  ( — ).  Calculated  at  x = 4.  48  In. , 
for  N = 81  (°),  for  N = 83  (x),  a = 0. 143  in. , p = 0.  25  in. 
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Figure  3-24.  Intensity  Behind  Hollow  Dielectric  Wedge,  k = 2.6,  With 
18  in.  x 18  in.  Walls,  Thicknesses  0.226  In.  and  0.232  in.  at  \=  1.259  in. 
Measured  with  dipole  probe  ( — ).  Calculated  at  x = 4.  48  in.  (x), 
at  x * 4. 61  In.  (°).  P = 0. 25  In. , a = 0. 136  in.  for  y < 0, 
a = 0. 130  in.  for  y > 0.  N = 81. 
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Figure  3-28a.  Computed  intensity  of  total  field  |Et|2  inside  hollow  90°  wedge  with 
1/8  in.  walls.  Dielectric  constant  2.6.  N=  161,  \=  1.259  in.  The 
arrows  show  spacings  of  nulls  predicted  bv  Equation  3-7. 

See  Figure  3-7  for  definition  of  o..  (Sheet  1 of  2) 
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Figure  3-29.  Intensity  and  phase  near  90 ; dielectric  wedge; 
sidewall  lengths  10  in.;  thickness  0.125  in. 
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1.  Intensity  Near  Hollow  Wedge  for  Figure  3-32.  Phase  Near  Hollow  Wedge  for 

22  Grey  levels;  0j  - 0°  22  Grey  levels;  0:  0® 


Figure  3-35.  Power  transmittance  of  dielectric  hemi -cylinder  for  A = 90°. 

Inner  diameter  4. 87  in. ; thickness  0. 125  in. ; wavelength 
1.259  in.  Moment  method  calculations.  Data  are  symmetric 
about  X = 0. 
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Figure  3-36.  Phase  delay  corresponding  to  Figure  3-35. 


Figure  3-37.  Power  transmittance  for  the  conditions  of  Figure  3-35,  except  with 
= 45°.  The  dashed  lines  show  the  co-ordinates  of  the  inner 
radome  surface. 
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Figure  3-40.  Power  transmittance  of  hemi-cylinder  with  dielectric  strips  for  = 90 : . 

Inner  diameter  of  hemi-cylinder  12.  38  cm;  thickness  0.  32  cm;  wave- 
length 3.  2 cm.  Computed  by  moment  method.  Radii  of  approximating 
cylinders  were  0.180  and  0.086  cm.  Data  are  symmetric  about  X = 0. 
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Figure  3-41.  |E  | Inside  Dielectric  Slab. 

at  y = 0,  for  N = 32,  a = 0. 071  in. , 


Thickness  0. 126  in. 
P=  0.126  in. 
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Figure  3-43.  |eT|  Inside  Dielectric  Slab,  Thickness  .235  in.  for  Grazing  Incidence. 

Parallel  Polarization 
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4. 1  APPARATUS  AND  PROCEDURE 

The  apparatus  is  sketched  in  Figure  4-1.  It  consisted  of  a signal  generator,  a network 
analyzer,  and  a probe  positioner.  The  probe  antenna  was  either  a half-wave  dipole  or 
an  open  end  waveguide.  The  transmitting  antenna  was  20'  from  the  probe  so  the  field 
was  approximately  that  of  a plane  wave.  Measurements  were  made  by  moving  the  probe 
on  the  paths  shown  in  Figure  2-2,  with  the  dielectric  scatterer  present  and  then  after 
the  scatterer  was  removed.  The  difference  between  measured  values  of  intensity  is 
the  power  transmittance,  which  is  normalized  to  the  free  space  value;  similarly  the 
difference  in  phase  values  is  the  argument  of  the  complex-valued  transmittance. 


4.2  RESULTS  FOR  PERPENDICULAR  POLARIZATION 
4.2.1  SLAB;  4 IN  WIDTH;  NORMAL  INCIDENCE 


For  the  slab  with  thickness  0.  24  in. , measured  values  of  I E /E  |“  are  shown  in 
Figures  3-6  and  3-7  for  two,  distinct  probes.  One  was  an  open-end  waveguide.  The 
other  was  a half-wave  dipole.  The  measured  data  depend  on  the  probe,  but  the 
differences  between  the  results  decrease  with  distance  from  the  slab.  For  either  probe, 
the  discrepancies  between  measured  and  computed  results  are  smaller  in  Figure  3-7 
than  in  Figure  3-6.  This  result  is  reasonable  since  N was  larger  for  Figure  3-7  than 
for  3-6.  In  Figure  3-7  discrepancies  are  smaller  for  the  dipole  probe.  Therefore  the 
dipole  seems  to  be  preferable  to  the  open-end  waveguide. 


Figure  3-8  shows 
0. 125  in. 


eVi2 


measured  with  a dipole  probe  for  a slab  with  thickness 


4.2.2  SLAB:  10  IN.  WIDE,  NORMAL  INCIDENCE 
i T 1 , 2 

Figure  3-10  shows  I E /E  I shows  measured  values  for  a slab  thickness  0.  25  in.  The 
dipole  probe  gives  the  variation  of  the  computed  curves,  and  discrepancies  are  smaller 
for  the  dipole  than  for  the  waveguide  except  in  the  center  of  the  graphs  for  Y = 0.  75  X. 


T 1 1 2 

Figure  3-11  shows  |E  /E  ' of  the  thinner  slab,  0.125  in.  thick,  for  a dipole  probe. 
Discrepancies  are  very  small.  The  reason  is  that  the  small  values  of  a gave  accurate 
computations.  In  contrast  for  Figure  3-10  N was  small,  and  a was  too  large  for 
accurate  computations.  This  approximation  is  further  discussed  in  Section  5. 
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4.  2.  3 SLAB  WITH  DIELECTRIC  STRIPS 


i T 1 1 2 

Measured  values  of  lE  ,/E  1“  are  shown  in  Figure  3-13.  The  discrepancies  increase 
with  increasing  X.  The  discrepancies  are  small  (less  than  0. 4 dB)  for  x > 0,  but  as 
much  as  0.  8 dB  for  x < 0.  This  asymmetry  results  because  the  computations  were 
done  with  the  strips  at  positions  asymmetric  about  the  slab  center.  In  Figure  3-13 
the  strips  near  the  fifteenth  more  nearly  represent  the  experimental  model  than  do 
those  near  the  first,  which  is  on  the  side  x < 0. 


4.2.4  SLAB,  45°  INCIDENCE 

i T 1 1 2 

Figure  3-17  shows  measured  values  of  intensity  I E /El  on  a transverse  path,  which 
had  a fuxed  x value  in  Figure  2-2b.  The  transition  between  the  directly  illuminated  and 
shadowed  regions  is  apparent.  The  discrepancies  are  large  in  the  shadowed  region 
because  the  cells  were  large;  with  a = 0. 143  in.  This  case  corresponds  to  that  in 
Figure  3-10;  only  the  incidence  angle  and  probing  paths  differ. 


4.  2.  5 SLAB,  GRAZING  INCIDENCE 

T 12 

For  grazing  incidence,  Figure  2-2c,  Figure  4-2  shows  Ie  /E  1“  measured  at  three 
values  of  the  spacing  between  the  probe  (a  half-wave  dipole)  and  the  slab.  The  spacings 
were  0. 1,  0.  2,  and  0.  3 in.  Figure  4-3  shows  |E  ^/E1!  ^ at  spacing  0.  4 in.  Notice 
that  spacing  is  less  than  the  y co-ordinate  by  half  the  slab  thickness.  For  example, 
spacing  0.  4 in.  corresponds  to  Y equal  0.  525  in.  Because  of  tolerances  in  slab  thick- 
ness the  value  of  y may  be  in  error  by  as  much  as  0.  01  in. 

4.  2.  6 HOLLOW  WEDGE,  0.  23  IN.  THICK,  AXIAL  INCIDENCE 

Figure  3-23  (and  3-24)  show  lEi/E  I for  a hollow,  90' wedge  formed  of  two  orthogonal 
slabs.  Wall  thicknesses  were  0.  226"  and  0.  232". 

4.3  RESULTS  FOR  PARALLEL  POLARIZATION 


The  probe  influences  measurements  for  parallel  polarization  when  it  is  near  a slab 
because  the  field  has  a considerable  evanescent  constituent.  That  is,  the  guided  wave 
amplitude  decreases  rapidly  with  increasing  distance  from  the  slab.  Since  the  probe 
spans  a finite  width,  it  averages  the  field.  Moreover,  the  response  to  the  field  depends 
on  probe  orientation  for  an  unbalanced  probe.  For  example,  Figure  4-4  shows 
measured  values  of  Ie^I^  for  two  orientations  of  the  dipole  probe  near  a slab,  thick- 
ness 0. 125  in.  The  orientations  are  approximately  symmetric,  with  the  probe  merely 
inverted,  yet  the  data  differ  considerably.  Figure  4-5  shows  results  for  a slab  0.  60  cm 
thick.  The  average  of  the  Ie^I  - values  for  the  two  orientations  is  shown  in  Figure  4-6. 
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The  data  in  Figure  4-6a  should  be  compared  with  those  in  Figure  3-42.  The  value  of 
Y we  expect  is  the  sum  of  half  the  probe  width  (0.  313  in. ),  half  the  slab  thickness 
(0.  063  in. ),  and  probe  clearance  from  the  slab  (0.  040).  The  sum  is  0.  416  in.  In 
Figure  3-42  the  graphs  1.  02  and  1.  27  cm  resemble  the  average  of  the  measured  data 
in  Figure  4-6. 

Figure  4-7  shows  deep  minima  near  a slab  0.432  in.  thick. 
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Figure  4-1-  Apparatus  for  Transmittance  Measurements.  NA 
represents  a network  analyzer,  G a generator,  T a transmitting 
antenna,  and  P a probe  antenna 
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Figure  4-2.  Intensity  Near  Slab  for  Grazing  Incidence. 
Dielectric  constant  2.  6;  thickness  0.  23  in. 

The  dashed  line  is  the  level  without  the  slab. 
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Figure  4-3.  Measured  Intensity  Near  Slab  at  Y = 0.  525  in. 
The  dashed  curve  is  the  level  without  slab  thickness  . 23  in. 
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Figure  4-7.  Measured  Power  Transmittance  Near 
0.  432  in.  Thick  Dielectric  Slab  for  Grazing  Incidence,  Polarization  Parallel 
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5.  DISCUSSION 


5.  1 PERPENDICULAR  POLARIZATION 

5. 1.  1 SLABS,  NORMAL  INCIDENCEj  INTERNAL  FIELDS  AND  AN  ACCURACY 
CRITERION 

T 

The  values  of  |E  I shown  in  Figures  3-4  and  3-5  oscillate  as  functions  of  X.  The 
values  on  the  slab-midplane  vary  about  the  value  for  an  infinitely  broad  sheet.  The 
exact,  analytical  value  for  the  mid-plane  is 

EP  =tQ1  (l-r10expkd’)  1(expikd'/2)  (11) 

where  t^  and  r^Q  are  Fresnel  interface  transmittance  and  reflectances,  and  d'  is 

dvK-sin“i,  where  d is  thickness,  and  i is  incidence  angle.  The  moment  method  results 
also  depend  on  y,  as  expected. 

In  the  following  discussion,  we  use  the  term  "discrepancy"  for  the  difference  between 
measured  and  computed  values. 

T 12 

Figure  3-6  shows  I E /E  I measured  once  with  an  open-end  waveguide  and  once  with 
a half-wave  dipole.  The  shapes  of  the  curves  for  the  dipole  more  closely  resemble  the 
computed  values  than  do  those  for  the  waveguide.  Discrepancies  are  also  smaller 
except  near  the  center  of  the  curve  for  Y = 0.  75  X.  The  discrepancies  are  large  for 
either  antenna. 

The  large  discrepancies  in  Figure  3-6  result  because  the  cell  size,  0.  114  X,  violates 
the  assumption  that  the  field  be  constant  in  each  cell.  For  a smaller  cell  size, 

0.  042  X,  Figure  3-7  shows  smaller  discrepancies.  In  this  figure,  discrepancies  are 
larger  for  Y = 0.  5X  than  for  Y = 1.  OX;  apparently,  the  probes  disturb  the  field  more 
when  they  are  nearer  the  slab.  Note  too  that  discrepancies  are  smaller  for  the  dipole 
probe.  Therefore  the  dipole  seems  better  for  perpendicular  polarization.  Apparently 
it  has  better  spatial  resolution  than  the  open-end  waveguide. 

Figure  3-8  shows  small  discrepancies  for  a thinner  slab.  In  this  case  calculations 
were  done  with  a equal  0.  715  in.  Again,  discrepancies  are  larger  for  the  probe 
nearer  the  slab. 
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Now  we  consider  a quantitative  accuracy  criterion  for  cell  radii.  We  plot  the  maximum 
discrepancy,  in  a graph  of  I E^/E*l“  for  each  y as  a function  of  optical  path  length 
corresponding  to  cell  diameter;  that  is, 


I 

0Q=  360o(2a  (12) 

Consider  the  data  measured  with  the  dipole  and  shown  in  Figures  3-6,  3-7,  and  3-8. 

The  values  of  for  the  respective  figures  are  132L,  49°,  and  66°.  Figure  5-1  shows 
maximum  discrepancies  plotted  as  functions  of  90  for  two  distances  between  probe 
and  slab  mid-plane. 

Figure  5-1  suggests  that  smaller  correspond  to  smaller  discrepancies.  For  the 
data  points  corresponding  to  the  two  smaller  values  of  0Q,  and  for  which  discrepancy 
magnitudes  are  relatively  small,  the  data  for  Y = X/2  seem  more  accurate  than  those 
for  Y = X . Apparently  the  probe  affects  the  measuremements  more  when  it  is  closer 
to  the  slab  than  it  does  at  larger  distances. 

We  conclude  that  maximum  discrepancies  of  reasonable  magnitude  (approximately  10CTc 
in  power)  can  be  obtained  with  0O  equal  66c.  This  value  results  with  a = 0.  042  X for 
k equal  2.6. 

5.1.2  SLAB:  25.  4 CM  WIDTH,  NORMAL  INCIDENCE 

For  the  thicker  slab.  Figure  3-10  shows  that  the  dipole  gives  somewhat  better  agree- 
ment with  the  shape  of  the  computed  curve  than  does  the  waveguide.  Discrepancies 
are  large  for  either  probe.  For  the  thinner  slab,  discrepancies  are  small  for  the 
dipole  probe;  see  Figure  3-11.  The  maximum  discrepancies  are  plotted  in  Figure  5-2 
as  functions  of  <pQ.  By  comparing  Figures  5-1  and  5-2,  we  see  discrepancies  are 
smaller  for  the  broader  slab.  If  we  choose  0.  6 dB  as  a criterion  then  for  q>Q  = 66°, 
this  discrepancy  can  be  achieved  at  Y = X/2  for  the  broader  slab;  however,  this  value 
requires  Y be  greater,  approximately  X. 

5. 1.  3 GUIDED  WAVES:  AN  INDEPENDENT  MODEL;  EXAMPLES  FOR  NORMAL 
INCIDENCE 

The  measured  and  computed  graphs  show  oscillations.  Examples  are  numerous.  The 
internal  field  magnitudes  in  Figures  3-4  and  3-5  vary  about  the  value  for  a sheet  of 
infinite  width.  At  normal  incidence,  Figure  3-11  shows  small,  but  distinct  oscillations, 
and  Figures  3-16  and  3-21  show  wide  variations.  These  oscillations  occur  for  both 
measured  and  computed  values.  Their  explanation  is  interference  of  guided  waves  and 
plane  waves  that  propagate  into  and  through  the  slabs  and  wedges.  This  interpretation 
follows  from  the  following  model. 
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For  notational  simplicity  let  us  start  with  external  fields  and  consider  a slab.  Suppose 
the  total  field  in  the  region  shadowed  by  a slab  consists  of  a plane  wave  and  a slab- 
guided  wave.  The  field  is 


T 

E = Ae 


-vd  ik  s 
e g 


^ ik  s cos  a 
T_,_  e o , 


(13) 


2 2 1/2 

where  s is  path  length  along  the  slab  or  (X  + Y“)  , kQ  is  2t  / \0,  with  \Q  the  free 

space  wavelength,  A is  the  guided  wave  amplitude  and  T±  the  complex -valued  plane 
wave  transmittance  of  a sheet  with  properties  of  the  slab.  The  symbol  d represents 
normal  distance  from  the  slab,  kg  is  the  propagation  constant  of  the  slab,  and  v is  the 
decay  factor.  The  angle  a is  the  angle  between  the  plane  of  the  slab  and  the  incident 
wave  normal.  The  solution  of  the  boundary  value  problem  for  slab  guided  waves  gives 
u and  v. 

From  symmetry  we  assume  the  guided  wave  to  have  even  variation  about  the  midplane. 

To  determine  k^  and  v we  must  solve 
© 

(ua/2)  tan(ua/2)  = (va/2)  (14) 

where  a is  slab  thickness,  and 


2 2 2 

u = k k - k , (15) 

o g 

with  the  condition 


2 2 2 

u +v  =k  {k-  1)  (16) 

o 

To  determine  u and  v,  graph  the  left  side  of  Equation  14  as  a function  of  ua/2  and  find 
the  intersection  with  the  circle  of  Equation  15,  when  written  as  (ua/2)"  + (va/2)  = 
kga^(K-l)/4.  For  example,  Figure  5-3  shows  a graph  for  2a  = 0. 125  in. ; \0  = 

1. 259  in. , and  k = 2. 6.  We  find  ua/2  is  0. 118  tt  so  kg  is 


k =1.73  /in.  (17) 

g 


and  v = 0.  72  /in. 


(IS) 
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As  a check  on  our  model,  we  shall  examine  spacings  of  fringes  as  in  Figure  3-16.  The 
condition  for  minima  spacings  is  that  the  two  terms  on  the  right  side  of  Eqatuion  13 
differ  by  2 t in  distance  s.  Thus 

6s=2ff/(k  -k  cos  a)  (19) 

g o 


To  derive  spacings  inside  the  slab  we  must  utilize  the  internal  refracted  and  reflected 
waves.  For  the  orientation  in  Figure  3-1,  assume  the  field  is 


_ i(k  x 
T x 

E = T e 
d 


+ V) 


i(k  x - k y)  ik  x 

Re*  y + Ae  g 
a 


(20) 


where  Td  and  Rd  are  transmittance  and  reflectance  for  the  internal  waves  (assumed 
plane),  =vxk0  sin  r,  ky  = \ Kk0  cos  r,  where  r is  the  angle  of  refraction.  Consider 
the  midpfane,  so  y is  zero.  Recall  Snell's  law,  that  sin  r = sin  i,  with  i the 
incidence  angle,  and  i + a = ir/2.  Thus,  Equation  20  gives 


ik  x cos  a ik  x 

ET  = (T  +R  )e  ° +Ae  8 (21) 

d d 

Equation  21  has  the  form  of  Equation  13  so  fringe  spacings  inside  and  outside  are 
equal.  This  result  is  physically  intuitive. 

The  model  assumed  for  Equation  13  was  one  plane  wave  and  one  guided;  but  this  model 
is  too  simple.  For  example,  consider  normal  incidence  on  a slab,  Figure  2-2a. 

The  oscillations  of  ET  inside  the  slab  can  be  explained  by  interference  of  guided  waves 
and  plane  waves,  but  from  symmetry  we  must  assume  two  guided  waves  that  travel  in 
opposite  directions.  We  also  assume  the  internal  field  contains  plane  waves,  which  for 
normal  incidence  are 

ik  Jk  y -ikQ\/ic  y 

E = T e ° + R ,e  , (22) 

p d d 

For  simplicity  let  y = 0 so 

E = T + R = S (23) 

p d d 
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Thus  the  internal  field  is 


I 


_ ik  x -ik  x 

ET  = S + Ae  S + Be  g 


(24) 


and  intensity  is 


Ie  I “ = C + 2 AB  cos  ^2k^x  + (a-/3)J  + 2|ASlcosjk^x  -1-  (a-cr)J 
+ 2 iBSl  cos  £k^x  + (cr-/3)J 

2 2 2 

where  C = I A I + B“  + |S|  , and  a,/3  and  <r  contain  phase  of  A,  B and  S. 
F rom  Equation  25  we  obtain  fringe  spacings. 


AS  = 7T/k  , 
g 


and 


(25) 


(26) 


AS  = 2 ff/k  . (27) 

g 

For  the  data  in  Figures  3-4  and  3-5,  we  have  kg  = 1.  96  ir/inch.  Therefore  AS  from 

Equation  27  is  1.  02  in.  This  value  agrees  closely  with  the  1.  00  in.  spacings  computed 

by  the  moment  method;  see  Figures  3-4  and  3-5.  This  spacing  was  determined  by 

averaging  the  spacings  in  the  graph.  Note  that  Equation  27  results  from  the  terms  with 

S as  a factor.  Physically,  this  spacing  corresponds  to  interference  of  plane  and  guided 

waves.  The  smaller  spacing,  Equation  26,  arises  from  interference  of  oppositely 

directed  waves.  This  smaller  period  is  apparently  absent.  The  reasons  are  that  A 

and  B are  smaller  than  S,  and  that  we  have  two  terms  with  period  2 7r/k  in  Equation  27. 

8 

The  mechanism  for  guided  wave  excitation  may  be  the  larger  phase  shift  over  the  slab 
and  negligible  phase  shift  away  from  the  slab.  The  wave  is  refracted  into  the  slab  near 
the  ends;  that  is,  the  wavefront  is  bent  into  the  slab  near  its  ends. 
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For  the  broader  slab  we  have  more  spacings  than  in  Figure  3—4.  Figure  3-9  shows 
spacings  from  moment  method  calculations  to  be  1.  0 in.  for  the  1/4  in.  thickness 
and  1.  14  in.  for  the  1/8  in.  thickness.  From  Equation  27  we  have  1.  02  in.  as  stated 
above  for  1/4  in.  thickness  Ov  = 1.  96  7r/inch)  and  1.  18  in.  for  1/8  in.  thickness 
(Iq,.  = 1.  70  7r/inch).  The  agreement  between  the  moment  method  and  Equation  27  is 
good. 

Now  consider  external  fields  in  Figures  3-10  and  3-11.  It  is  clear  that  the  moment 
method  agrees  with  measurement  for  the  thinner  slab  (Figure  3-11);  for  Figure  3-10 
discrepancies  are  large  because  a was  large. 

Equation  27  also  gives  predictions  that  agree  with  measurement.  For  the  1/4  in. 
thickness,  Equation  27  gives  AS  as  1.  02  in.  ; the  average  spacing  of  8 measured 
fringes  in  Figure  3-10  is  1.  0 in.  For  the  1/8  in.  thickness.  Equation  27  gives  AS 
as  1.  18  in.  ; the  average  spacing  of  7 measured  fringes  in  Figure  3-11  is  1.  14  in. 

It  is  significant  that  fringe  spacings  depend  on  slab  thickness.  We  can  imagine  a simple 
model  in  which  the  slab  edges  scattered  free  space  type  rays.  In  this  case,  no  thick- 
ness dependence  would  exist.  Moreover,  the  fringe  spacing  would  be  given  by  the 
intensity  arising  from  an  expression  such  as 

ik  x -ik  x ik  v 

E = Ae  ° + Be  ° * Te  ° (28) 

where  for  emphasis  kQ  is  the  free  space  propagation  constant.  From  Equation  28  we 
obtain  periods 

AS  = */k  (29) 

o ' 

and 


AS  = 2ir/k  (30) 

o 

From  Equation  30  the  spacing  is  \ or  1.  259  in.  This  value  exceeds  the  1.  14  in.  value 
for  1/4  in.  thickness  and  the  1.  0 in.  value  for  the  1/8  in.  thickness. 

We  conclude  edge  diffracted  rays  are  absent,  but  that  guided  waves  axist  and  they  might 
be  considered  to  arise  from  refraction  at  the  edges.  We  would  call  these  guided  rays 
in  a geometric  diffraction  theory  context.  Guided  waves  also  correspond  to  the  rapid 
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decrease  of  the  amplitude  of  the  oscillations  with  distance  in  Figures  3-10  and  3-11; 

these  slabs  were  relatively  broad  with  widths  8 wavelengths.  For  the  narrower  slabs, 

I T T,  9 

the  rapid  oscillations  of  IE  /El“  near  the  edges  may  suggest  cylindrical  waves  from 
the  edges,  but  we  have  not  examined  this  possibility. 

We  emphasize  that  the  model  leading  to  Equation  13  is  independent  of  the  moment 
method  yet  it  yields  spacings  that  agree  with  those  of  the  moment  method.  In  addition, 
both  the  moment  method  results  and  Equation  13  yield  results  that  agree  with  measured 
data.  The  evidence  for  guided  waves  is  strong;  we  have  identified  a simple  wave  type , 
guided  waves,  in  the  numerical  results  from  the  moment  method. 

3.  1.4  SLAB  WITH  DIELECTRIC  STRIPS 

The  calculated  values  in  Figure  3-14  show  the  strips  increased  transmittance  averaged 
over  the  slab  by  approximately  1/3  dB.  In  Figure  3-13,  discrepancies  between 
measured  and  computed  values  are  small  for  X > 0,  but  larger  for  negative  X.  These 
discrepancies  arise  because  the  experimental  model  differed  from  the  computational. 
The  strips  locations  for  X > 0 agree  for  both  models.  In  Figure  3-12  the  configuration 
near  the  15th  strip  is  that  of  the  experiment  model;  however,  near  the  1st  the  compu- 
tational model  differs  from  the  experimental. 

5.1.5  SLAB.  45  DEGREE  INCIDENCE 

The  deep  minima  in  Figure  3-15  can  be  explained  by  interference  between  guided  waves 
and  the  internal  plane  waves.  From  Equation  19,  with  kg  equal  1.  96  tf/inch  and  a 
equal  45°  we  obtain  AS  = 2.  39  in.  which  agrees  well  with  the  2.  4 in.  spacing  in 
Figure  3-19.  The  shallower,  more  rapid  oscillations  apparently  result  from  inter- 
ference of  the  guided  wave  travelling  in  the  positive  X direction  of  Figure  3-1  and 
a smaller  amplitude  guided  wave.  The  second  term  in  Equation  25  represents  such 
waves,  and  Equation  26  gives  fringe  spacing  0.51  in.  This  value  agrees  with  the  period 
of  the  rapid  oscillations.  Recall  that  cells  are  spaced  by  0.  25  in.  and  in  Figure  3-15 
the  period  corresponds  to  two  cells.  The  out  of  phase  behavior  of  the  small  oscillation 
produced  by  adding  a cell  also  supports  the  interpretation  of  guided  wave  interference. 
The  added  phase  on  reflection  is  2 k„.  P or  2(1.  96;r)/4  = 0.  dSn.  The  phase  is  almost 
reversed  by  adding  a cell. 

Notice  that  the  fringe  pattern  is  established  even  for  N = 10  which  corresponds  to  slab 
length  only  2 wavelengths. 

The  external  field  at  distance  0. 1 in.  from  the  slab  shows  fringes  with  spacings  closely 
approximating  the  2.4  in.  value  predicted  by  Equation  19. 
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For  the  transverse  path  of  Figure  2- 2b,  (Figure  3-17),  the  abrupt  decrease  for  Y 
approximately  4 in.  fits  the  notion  of  a guided  wave. 

5. 1.  6 FLAT  SLAB,  GRAZING  INCIDENCE 

Figure  3-19  shows  that  the  computed  internal  field  intensity  exceeds  the  incident  field 
intensity.  The  field  magnitude  I E X I reaches  values  almost  twice  that  of  |e*|.  Again 
we  see  the  out-of-phase  behavior  produced  by  adding  one  cell.  Figures  3-21  and  3-22 
show  that  the  field  changes  rapidly  with  perpendicular  distance.  By  comparing  Figures 
3-20  and  3-21  with  Figures  4-2  and  4-3  we  see  excellent  agreement  between  measured 
results  and  those  from  the  moment  method.  Notice  in  Figure  3-20  that  the  small  rapid 
oscillations  have  period  0. 5 in.  in  the  region  of  the  slab,  0 in.  s x s 10  in. ; however, 
for  x > 0 the  period  exceeds  0. 6 in.  Thus  we  have  a free  space  reflected  wave  for 
x < 0 and  interfering  guided  waves  for  X > 0. 

5.1.7  HOLLOW  WEDGE 

The  internal  field,  shown  in  Figure  3-22,  for  thickness  0.  25  in.  at  axial  incidence  has 
oscillations  like  those  of  the  slab,  Figure  3-15.  Again,  we  have  interference  between 
guided  and  free  space  plane  waves.  In  Figure  3-22,  the  amplitude  oscillates  between 
0. 1 and  1. 1,  about  the  value  0.  6 for  the  internal  plane  waves.  Therefore,  the  guided 
wave  amplitude  is  estimated  at  0.5.  This  estimate  also  applies  for  the  slab  at  45° 
incidence;  see  Figure  3-15. 

Figures  3-23  and  3-24  show  that  distinct  thicknesses  for  the  two  walls  of  the  wedge 
introduce  small  fluctuations  into  the  field  on  the  transverse  plane.  These  fluctuations 
appear  in  the  measured  data. 

The  intensity  I E /E  I on  the  symmetry  plane,  Figure  3-25,  is  influenced  by  waves 
that  are  scattered  by  the  edges  of  the  wedge.  The  edges  at  the  far  ends  of  the  walls 
project  to  x = 7. 1 in.  in  the  co-ordinate  system  of  Figure  2- 2d,  and  a third  edge  is  at 
the  vertex.  Assume  three  rays  as  in  Figure  5-4.  Let  the  fields  of  these  three  rays 
be  a sum  of  three  plane  waves 


The  intensity  is 


I 

I 


2 9 2 r 

I E I = |T  f - a + 4|aT  | cos  (k-k  ) x + t - 
e p1  p L x 


where  t is  the  phase  of  T ando;  that  of  a.  From  Equation 


(5x)e  = (1  - cos 0)  1 


«]  (33) 

33  the  period  of  |e  1“  is 
(34) 


which  is  plotted  in  Figure  5-5.  For  0=  90°,  ( x)g  is  X.  In  Figure  3-25  the  spacings 
for  5 in.  Sin.  are  approximately  X . Furthermore,  fringe  spacings  increase 

with  x,  this  trend  is  predicted  by  Equation  34,  as  shown  in  Figure  5-5.  The  variation 
for  x < 4 in.  are  small  suggesting  the  edge  diffracted  waves  vary  with  position. 

The  field  structure  behind  a wedge  is  evidently  complicated,  for  it  includes  guided 
waves,  edge- diffracted  waves,  and  waves  propagated  through  the  walls.  Figure  3-27 
suggests  that  the  field  is  approximately  a plane  wave  because  the  equiphase  surfaces 
are  relatively  flat.  However,  Figure  3-26  shows  substantial  intensity  variation 
especially  near  the  wedge  walls. 

Equation  19  predicts  the  fringe  spacings  for  the  hollow  wedge  with  0. 125  in.  thick 
walls.  Table  5-1  summarizes  the  data.  For  the  wedge,  we  have  two  values  of  incidence 
angle  for  each  0^  because  we  have  two  slabs  forming  the  wedge;  of  course,  for  0j  = 0°, 
the  two  values  of  a are  equal,  45°. 


Table  5-1 

Fringe  Spacings  From  Equation  7 for  k^  = 1.  73  7r/inch,  2a  = 0. 125,  k = 2.  6,  A = 1.  259 
a 


0i 

(deg) 

a 

upper  slab 
(deg) 

O! 

lower  slab 
(deg) 

<5s 

upper  slab 
(inch) 

<5s  j 

lower  slab 
(inch) 

22.5 

22.5 

67.5 

7.407 

1.  7S4 

45 

0 

90 

14.29 

1.156  | 

67.  5 

22.5 

67.5 

7.407 

1. 784  | 

0 

— 

45 

45 



3.30 

3.30 

Figures  3- 28a  and  b have  these  spacings  superimposed  on  the  moment  method  calcu- 
lations of  | ETI 2 . We  observe  agreement  is  excellent.  We  also  see  the  rapid 
oscillations  predicted  by  Equation  24. 
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To  study  the  exponential  decrease  of  guided  waves  we  computed  |E  J.  Figure  5-6 
shows  ES  on  two  Lines,  which  intersect  the  x axis  at  4 in.  and  at  5 in.  The  co- 
ordinate is  parallel  to  the  incident  wavefront.  In  this  case,  = 45°.  The  graphs 
show  large  values  for  greater  than  2 in.  or  3 in. ; that  is,  near  the  upper  wall  of 
the  wedge.  The  figure  also  shows  Ae~vc^,  where  A is  the  value  from  the  moment 
method  calculation  just  outside  the  wall.  The  exponential  decrease  agrees  fairly  well 
with  the  moment  method  results.  The  approximately  constant  values  of  lE^I  for  Y less 
than  approximately  2 in.  arise  mainly  because  of  the  lower  half  of  the  wedge. 

S T 

Although  the  scattered  field  E is  theoretically  interesting,  the  total  field  E is 

directly  observable  and  of  interest  for  radomes.  Therefore,  we  show  |E^|  2 for  the 

scattered  field  of  Figure  5-6  in  Figure  5-7  and  the  phase  or  arg  E^  in  Figure  5-8. 

In  Figures  5-7  and  5-8  the  dashed  lines  show  the  values  for  an  infinitely  broad  sheet. 

Near  the  upper  slab,  the  moment  method  results  differ  very  much  from  those  for  the 

flat  sheet.  Away  from  the  upper  slab  and  behind  the  slab  illuminated  at  normal 

incidence,  the  moment  method  results  oscillate  about  the  flat  sheet  results.  The 

oscillations  imply  the  existence  of  a field  distinct  from  that  for  a flat  sheet.  In  such  a 

description,  we  would  exclude  the  guided  wave  as  a constituent  combining  with  the 

wave  propagated  through  the  lower  slab  because  the  exponential  decay  is  strong. 

Instead,  we  might  expect  scattering  from  the  tip  or  the  far  ends  of  the  wedge.  Of 

course  the  guided  waves  may  be  scattered  by  the  edges  or  ends  of  the  dielectric 

structure. 

g 

We  obtain  a comprehensive  view  of  the  field  near  the  slab,  from  the  displays  of  E 
by  a computer-driven,  oscilloscope  plotter.  Photographs  of  I ES I “ and  the  phase  of 
E®  are  shown  in  Figures  3-29  and  3-30.  These  figures  show  the  highly  excited  guided 
waves  for  grazing  or  nearly  grazing  incidence,  and  the  relatively  large  values  of  os 
for  such  incidence.  The  figures  show  the  smaller  value  of  <5s  for  0i  = 45°.  The  less 
visible  fringes  apparently  occur  because  of  interference  between  waves  from  the 
vertex  and  slab  ends  with  the  plane  wave  through  the  slabs.  The  wave  processes 
producing  the  weaker  fringes  may  be  significant  for  boresight  error  with  an  aperture 
antenna  because  many  of  these  antennas  measure  fields  through  two  balanced  halves 
of  the  antenna.  Small  asymmetries  can  unbalance  the  antenna. 

5.1.8  HOLLOW  WEDGE:  BORESIGHT  ERROR 

Although  few  radomes  have  the  shape  of  a wedge,  it  is  useful  to  estimate  the  boresight 
error  of  this  idealized  shape.  The  method  is  to  determine  the  tilt  of  the  wavefront, 
as  described  in  Appendix  A.  To  compute  tilt,  we  fit  a least  squares  line  to  phase  delay 
data  like  that  in  Figure  5-8.  In  particular,  we  choose  a line  parallel  to  the  incident 
wavefront  for  Pj  = 0,  22.5°,  45°,  and  67.5°;  this  line  intersects  the  symmetry  axis  at 
x = 4 in.  The  line  represents  the  position  of  an  antenna  aperture.  For  each  angle  we 
sample  the  phase  at  a set  of  points  spaced  by  X0/2  and  fit  a least  squares  line  to  these 
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phase  values.  To  estimate  effects  of  aperture  size  we  use  nine  points  to  represent  a 
four  \Q  wide  aperture  and  seven  points  to  represent  a three  XQ  wide  aperture.  Typical 
results  for  = 22.5°  are  tilt  of  9.  2°A0  or  1.5°  boresight  error,  for  the  4AQ  aperture. 

Boresight  error  is  summarized  in  Figure  5-9.  The  graphs  have  the  shape  of  typical 
boresight  error  curves  with  a maximum  at  angles  such  that  the  projection  of  the 
aperture  begins  to  clear  the  tip  of  the  radome.  The  boresight  error  is  larger  for  the 
larger  aperture.  This  result  is  reasonable  when  we  recall  the  abrupt  change  in  phase 
near  the  wall;  see  Figure  5-8.  A smaller  aperture  misses  much  of  the  phase  delay 
increase.  Figure  5-9  also  shows  boresight  error  when  the  phase  is  computed  by  flat 
sheet  theory.  In  Figure  5-8,  these  values  are  the  straight,  horizontal  line  with  zero 
tilt  so  boresight  error  is  zero  for  = 45°. 

The  boresight  error  magnitudes  for  the  phase  delays  computed  by  the  moment  method 
are  larger  than  those  based  on  flat  sheet  theory.  The  error  is  more  convergent  for 
the  moment  method,  where  convergent  means  that  the  normal  to  the  tilted  phase  line 
shifts  toward  the  radome  tip,  away  from  the  antenna  normal.  This  sense  of  shift  is 
reasonable  because  the  guided  wave  is  slow,  increasing  the  phase  delays  near  the  wall. 
This  effect  is  observed  in  axially- symmetric  radomes,  as  we  show  in  the  next  section. 

5. 1.  9 AXIALLY-SYMMETRIC,  THIN  SHELL 


The  moment  method  results  computed  for  the  wedge  resemble  those  observed  in  axially 
symmetric  shells.  For  example,  Figure  5-10  shows  power  transmittance  |Tl“ 
measured  in  the  region  bounded  by  an  axially- symmetric  shell  with  thickness  0.  029  A0  , 
dielectric  constant  3.2,  length  l2.8An,  and  base  diameter  6.7\n.  Gimbal  angle  was 

. | O U 

10°.  | T I is  shown  for  both  horizontal  and  vertical  polarization,  and  for  the  probe 
scanned  over  paths  at  three  distances  from  the  tip.  The  data  for  vertical  polarization 
show  |T 1 2 exceeding  0 dB  for  y less  than  zero;  this  is  the  side  nearer  the  radome 
wall  as  sketched  in  Figure  5-11.  The  transmittance  greater  than  0 dB  near  the  wall 
suggests  guided  waves. 


Boresight  error  of  the  axially- symmetric  shell  was  measured  with  an  aperture 
antenna,  diameter  4. 46  AQ;  see  Figure  5-12.  The  error  was  computed  with  a surface 
integration  method  in  which  the  integration  extends  over  the  first  Fresnel  zone  of  the 
incident  wavefront  at  the  radome  tip.  This  method  has  given  accurate  results  for 
half-wave  radomes,  but  for  the  thin  wall  radome  discrepancies  are  large.  The 
measured  curves  are  more  convergent  than  the  calculated,  where  convergent  means 
that  the  radome  tilts  the  tracking  null  toward  the  radome  axis,  away  from  the  bore- 
sight axis  without  radome.  In  terms  of  wavefronts,  convergent  implies  greater 
phase  delay  along  the  lower  \vall  for  the  situation  of  Figure  5-11.  This  result  agrees 
qualitatively  with  the  effect  expected  of  a guided  wave,  which  is  a slow  wave.  This 
increased  phase  delay  is  evident  for  a wedge  as  in  Figure  5-8  for  > 0. 


5-11 


Figure  5-12  shows  phase  measured  in  the  thin  axially  symmetric  shell.  It  also  shows 
data  computed  with  the  surface  integration  method.  The  data  are  given  on  two  lines. 
One  the  horizontal  diameter;  the  other  0.5  in.  above.  The  discrepancies  are  large 
for  either  polarization  although  the  same  surface  integration  method  has  given 
accurate  results  when  applied  to  a half-wave  radome  of  the  same  size  and  shape. 

For  horizontal  polarization,  the  tilts  of  the  measured  and  computed  equiphase  sur- 
faces differ  appreciably  in  Figure  5-13.*  The  large  difference  in  tilts  seems  to 
explain  the  large  discrepancies  for  the  E-plane  boresight  error  in  Figure  5-12.  The 
measured  phase  delays  are  generally  larger  than  the  computer-especially  near  the 
wall.  This  behavior  suggests  guided  waves. 

5. 1. 10  HALF  WAVE,  AXIALLY  SYMMETRIC  RADOME 

The  intensity  variations  in  Figure  3-25  resemble  variations  that  were, observed  in  an 
axially  symmetric  radome.  This  radome  was  a half-wave  thick,  highly  streamlined, 
and  small,  with  base  diameter  5.  2A0  and  length  10.  3X0. 

The  measured  intensity  in  the  radome- bounded  region  is  shown  in  Figure  5-14  for 
axial  incidence.  The  measurements  were  made  with  an  open-end  waveguide  probe  in 
an  interferometer  setup  which  had  a vertical  propagation  path.  The  radome  was  sup- 
ported by  three  thin,  wooden  dowels  in  order  to  reduce  scattering  by  supports. 

The  intensity  variations  for  the  smaller  values  of  x have  period  between  0.  6A0  and 
0.S>o.  This  value  compares  reasonably  with  (fix)g  = 0.6\o  for  0 = 130°. 

The  intensity  near  the  tip  exceeds  the  incident  value  and  varies  with  a period  approxi- 
mately a wavelength.  It  seems  plausible  that  the  high  intensity  results  because 
guided  waves  are  excited.  The  guided  wave  interpretation  is  supported  by  the  spacings 
of  minima  near  the  tip.  The  spacing  for  an  assumed,  slab-guided  wave  is 


(fix)  = 27r/[k  (cos0)  1 -k  ] 
g g o 


(35) 


where  9 is  the  angle  between  the  slab  and  the  free  space  wave  direction.  We  obtain 
(fixJo-  = 1. 1 in.  which  compares  well  with  the  1.  2 in.  ±0. 1 in.  spacings  in  Figure  5-13. 

O 

5.1.11  HEMICYLINDER 

T 2 

Figure  3-35  shows  that  the  moment  method  values  for  I E | vary  about  the  results 
by  approximating  the  radome  as  locally  flat. 

* 

These  graphs  have  been  reversed  from  those  in  Figure  5-10  because  the  gimbal 
angle  for  measurement  is  the  negative  of  that  for  computation. 
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By  comparing  Figures  3-35  and  3-40,  we  see  that  the  addition  of  anisotripic  gratings 
increases  power  transmittance  by  0.6  dB  on  the  average  over  the  diameter. 

The  large  variations  near  the  wall  suggest  guided  waves. 

5.2  PARALLEL  POLARIZATION 

The  moment  method  results  agree  well  with  measured  data,  but  the  measurements  with 
the  dipole  probe  are  averaged  values  (over  the  probe  length).  The  accuracy  of  fields 
measurements  seems  rather  low,  but  the  measured  data  verify  the  moment  method 
calculations.  The  lack  of  accuracy  is  in  the  position  of  the  measurement  when  the 
probe  is  near  the  wall.  The  computed  field  decreases  rapidly  with  distance  from  the 
wall;  the  probe  averages  this  evanescent  field.  The  computed  field  is  defined  mathe- 
matically at  a point. 

Figure  4-7  shows  deep  minima  exist  near  a slab  for  parallel  polarization,  just  as  they 
do  for  perpendicular.  However,  the  two  polarizations  require  distinct  thickness  for 
equal  fringe  spacings. 
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Figure  5-6.  Computed  values  of  scattered  field  magnitude  for  90°  hollow  wedge 
thickness  0. 125  in. , = 1.  259  in. , N = 161,  = 45°.  The  co-ordinate  YA  is 

on  a line  45°  to  the  Y axis;  it  is  parallel  to  the  incident  wavefront  for  the  upper  graph 
Ya  = 0 occurs  at  x = 4 in. ; for  the  lower  x = 5 in.  From  moment  method:  (-);  with 
exp(-vd):  (--). 
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Figure  5-8.  Phase  for  Figure  5-7 
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Figure  5-12.  Measured  and  computed  boresight  error 
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Figure  5-14.  Power  transmittance  of  ogival  radome;  inner  surface  of  tip 
at  x = 0.  Dielectric  constant  5.  5,  thickness  0.  222\;  with 
0.  026\  thick  coating  of  dielectric  constant  3. 0 
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APPENDIX  A 

WAVEFRONT  TILT  AND  BORESIGHT  ERROR 


This  appendix  shows  that  radome  boresight  error  is  related  to  the  tilt  of  the  wavefront 
in  the  radome-bounded  region. 

For  reception,  the  received  complex  amplitude  of  a rectangular  component  is 


,R 


-// 


FTE1  dA 


(A-l) 


where  F is  the  aperture  distribution,  E1  is  the  amplitude  of  an  incident  plane  wave 
amplitude,  T is  the  complex-valued  transmittance.  The  integration  extends  over  the 
region  such  that  F has  appreciable  magnitude,  and  T depends  on  the  angle  of  arrival 
of  the  incident  plane  wave  as  well  as  position  inside  the  radome. 

To  simplify  the  notation  let  us  consider  a symmetric  situation.  Assume  an  axially- 
symmetric  radome  with  its  axis  horizontal  and  the  antenna  normal  rotated  in  the  hori- 
zontal plane  as  in  Figure  A-l.  Establish  a rectangular  co-ordinate  system,  co- 
ordinates x,  y,  z,  attached  to  the  antenna,  and  let  the  antenna  be  in  the  x = 0 plane. 

Let  the  incident  wave  normal  be  in  the  x-y  plane.  To  maintain  symmetry  about  z = 0, 
consider  the  field  polarized  horizontally  or  vertically.  The  incident  wave  is  exp  ikyy 
where  kv  is  (2ir/\0)  sin(6  0).  with  63  the  off  axis  angle  of  the  wave.  Let  the  aperture 
distribution  be  uniform  in  its  finite  domain.  Let  T be  exp  iA<$,  where  A$  is  phase 
delay  so  the  amplitude  variation  is  ignored.  Introduce  a phase  taper  through  exp  iav 
to  describe  conical  scanning.  From  Equation  A-l  we  have 


E^  = JJ exp  i [ A$  + (k  + a)  y]  dA 


(A-2) 


Now  approximate  A$  with  a polynomial  so  A$  is  p(y,  z)  J-  r(v,  z)  with  p the  polynomial 
and  r a remainder.  If  we  collect  ky  into  the  linear  term  of  the  polynomial  we  have 

<t>  = A$  + k^y 


2 3 2 2 

= a + (b*  + k ) y + cy  + dv  + ez  y + fz 


(A- 3) 


if  we  expand  to  the  third  degree.  Note  that  odd  powers  of  z are  absent  because  of  the 
assumed  vertical  symmetry.  By  combining  Equations  A-2  and  A-3  we  have 


For  conical  scanning  systems  the  boresight  condition  is 


|ER  (Qf)| 2 = |ER(-a)|2 


(A-5) 


Thus  adjustment  of  ky  gives  equal  power  in  two  lobes.  A different  condition  holds  for 
monopulse  type  antennas.  This  adjustment  reduces  phase  variation  so  we  can  expand 
the  exponential  in  <$.  Similarly  a is  usually  small  so  its  exponential  can  be  expanded. 
The  remainder  term  also  can  be  expanded  if  it  is  small.  After  some  algebra  we 
obtain  from  the  condition 

6|E|2  = |ER  (a)|2  - |ER  (-or)l2  = 0 (A-6) 


that  k 

y 


at  boresight  is 


k 

yo 


- b'  - - (e  + 3d) 


(A-7) 


Boresight  error  is  the  angle  60  such  that  kvo  is  (2ir/\0)  sin(60)o,  or  for  small  angles 
(2  7r  /A 0)  (60  )G.  Thus  boresight  error 

(6  0)q  = - jb'  + i (e  + 3d)  R2j/(2  7r/x0)  . (A-8) 

From  Equation  A-8,  for  small  higher  order  terms,  in  e and  d,  boresight  error  is 
approximately 

(6  0)q  = -b'/(2rr/\o)  . (A-9) 

Thus  boresight  error  is  approximately  proportional  to  the  linear  phase  tilt,  where  tilt 
means  the  magnitude  of  the  expansion  coefficient  b'  in  the  polynomial  fitted  to  A$. 

In  practice  we  fit  a least  square  surface  to  A$.  Note  that,  although  many  approxi- 
mations were  made,  none  restrict  A$  from  realistic  values. 
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Figure  A-l.  Co-ordinate  system 
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